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Abstract. The three-body Schrodinger operator in the space of square 
integrable functions is found to be a certain extension of operators which 
generate the exponential unitary group containing a subgroup with nilpo- 
tent Lie algebra of length k+1, k = 0,1,. . . As a result, the solutions to the 
three-body Schrodinger equation with decaying potentials are shown to 
exist in the commutator subalgebras. For the Coulomb three-body sys- 
tem, it turns out that the task is to solve - in these subalgebras - the radial 
Schrodinger equation in three dimensions with the inverse power poten- 
tial of the form r'"'^. As an application to Coulombic system, analytic 
solutions for some lower bound states are presented. Under conditions 
pertinent to the three-unit-charge system, obtained solutions, with /c = 0, 
are reduced to the well-known eigenvalues of bound states at threshold. 

PACS: 03.65.Ge, 03.65.Db, 03.65.Fd 



1. Introduction 

The goal of the present paper is to demonstrate an analytical approach 
for solving the three-body Schrodinger equation with translation invariant 
decaying potentials. The three-body Schrodinger operator H (the Hamil- 
tonian operator, henceforth) is represented by the closure of operator sum 
T + V. The kinetic energy operator T is defined so that its closure, denoted 
by T as well, is a self-adjoint operator on the domain D{T) and acting in 
L2(/?9) by 



T = - Y,(2m,)-'Ai. (1.1) 

l</<3 

The constants (referred to as masses) m,- > (/ = 1,2,3), the Laplacian 
A,- is in three-dimensional vectors r,- = (xi,yi,Zi) e R^, with absolute value 
e [0, oo). The potential energy operator V is a scalar translation invari- 
ant operator of multiplication by V, where real function V fulfills several 
assumptions: (Al) 

1 
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V= J] Vijin-rj) and Vij^O as |r,-ry|^oo. (1.2) 

l<;<y<3 

(A2) Vij is of the differentiability class C^iR^) and it is analytic everywhere 
except, possibly, at = rj for i j. (A3) The operator V is assumed to be 
a symmetric T-bounded operator in the sense of Kato [KatSl] (see also 
[SimOO] and the citation therein), with its domain satisfying D(V) D D(T). 
This assumption ensures the self-adjointness of H on D(T). 

Let Vij be the gradient in vectors r, ^ = r, - rj with the absolute value r, j. If 
one defines the sum Yii<j^ii by G (Lemma 2), with its z-component, and 
Vn by G'lV (n = 0, 1, . . .), then there exists a subset 1^ c (§2.2, eq. (2.3)) 
such that operators G^, Vq = V,V\, Vk form a (/c- + 2)-dimensional nilpo- 
tent Lie algebra J{. in /^^ (see Theorem 1), whereas V^+p is the operator of 
multiplication by zero for all integers p = 1, 2, . . . For a particular Coulomb 
three-body system, the latter leads to a well-known observation that bound 
states exist whenever one of the three charges has a different sign (see, for 
example, [BD87, FB92, MRW92]), though this is not a sufficient condition 
of boundedness (see also §3.1 or, in particular, eq. (3.1)). 

After establishing the nilpotency of and thereby the existence of its 
commutator subalgebras c ^ (c = 0, 1 , . . . , /f -i- 1) we come to an impor- 
tant conclusion that the eigenvalue E of fl^ in I?{R^) is the sum of E^, the 
eigenvalues of H,^ in L^ih), where /c runs, in general, from to oo (eq. (2.5)), 
whereas the eigenvalue equation for can be decomposed, to some extent, 
into two equations in distinct subspaces so that their common solutions with 
respect to the eigenvalues were equal to E*^, eq. (3.5). Then the eigenvalue 
E^ equals E*^ plus the correction due to the Hughes-Eckart term, eq. (3.3). 
In case that V represents the Coulomb potential, these two equations are 
nothing but the separated radial Schrodinger equations in three dimensions 
with the inverse power potential of r~'^~^ type. Namely, one of our main 
goals is to demonstrate that solutions to the three-body Schrodinger equa- 
tion in a Coulomb potential are approximated by solving a one -dimensional 
second order ordinary differential equation in L^(0,oo;dr), 

u';i(r) + iB,-m+l)/r^]-[AJr'^^])u,,i(r) = 0, 1 = 0,1,..., (1.3) 

with some real constants and Bk, the latter being proportional to E^. The 
subspace 1^ yields Vp = for all integers p > 1 , as a: = 0; in particular, Vi = 
is a familiar relation known as the Wannier saddle point (see eg [SG87]). 
Under appropriate boundary conditions functions UQj(r) are expressed in 
terms of the confluent Whittaker functions, and the associated eigenvalues 
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are proportional to oc -Aq/(4r^), where integers n>l+l. These eigen- 
values represent the energies of bound states at threshold, §3.2. Based on 
perturbative arguments one deduces that the Hughes-Eckart term does not 
influence the ground state of H^, Proposition 1; see also [SimVO]. 

Unlike the previous case, admitting /c = 1 we are dealing with operators 
in the subspace Jl^ = [SI, SI], for V\ = Air~^ is the transition potential in the 
sense of [FL71] while Vp = for all integers p>2;Vii^0 ensures a higher 
accuracy in obtaining the eigenvalues + E^ from a series expansion. Ac- 
cording to Case [Cas50], for potentials as singular as or greater, there 
exists a phase factor - proportional to the cut-off" radius ro - that describes 
the breakdown of the power law at small distances r. Hence E^ = E^^{rQ) 
for /<■ = 1, 2, . . . On the other hand, is proportional to (-1)'^ (Corollary 4), 
which separates attractive potentials from the repulsive ones thus bringing 
in different characteristic aspects of eigenstates. 

The task to solve the radial Schrodinger equation entailing singular po- 
tentials has been of a particular interest for years ([Cas50, FL71, Spe64, 
Yaf74]) and still it draws the attention of many authors ([Gao08, GOK"'"01, 
ISl 1, MEFOl, RobOO]). If following [Cas50] or [Spe64], there is no ground 
state for these singular potentials as well as there are only a finite number of 
bound states, as abstracted eg in [Gao98, Gao99a, Gao99b]. So one should 
expect a more valuable contribution to eigenvalue expansion at higher en- 
ergy levels; for more details, see [RS78, §XIII]. Even so, several attempts 
have been initiated to find the ground state energy for a particular class of 
singular potentials (see eg [BBC80, NCU94]). 

The calculation problems of singular potentials are beyond the core of the 
present paper, though the results on the subject are highly appreciated for 
several reasons (see §3). We shall not concern ourselves with accomplishing 
the task to find general solution to eq. (1.3) but rather demonstrate a method 
for obtaining it in the Coulomb case; only the eigenvalues Eq as well as E\, 
in some respects, will be discerned for some illustrative purposes (§3.2.1- 
3.2.3). 



2. SlMILARTTY FOR THE THREE-BODY HaMILTONIAN OPERATOR 

Throughout the whole exposition, we shall exploit several Hilbert spaces 
of square integrable functions. The typical of them are: L^(R^) (as the base 
space), L^(R^) (as a subspace of translation invariant functions), L?-(Ik) (as 
a space over vector space 1^ c R^), and (so that H^'SiH^ = The 
norm and the scalar (inner) product in a given Hilbert space will be denoted 
by II • II and (•,•); whenever necessary, the subscripts identifying the space 
will be written as well. 
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2.1. Translation invariance. This section summarizes some requisite re- 
sults obtained from the translation invariance of the potential energy opera- 
tor V. 

We say e(t) is a representation of the group of translations in R , denoted 
E(3), in the space of (R^yiunctions f(ri) if it fulfills e(t)f(ri) = f(ri + 1), 
where vector / = (txJyJz)- A Taylor series expansion of /(r, -i- f) yields 
e{t)f{ri) = exp(if • Pi)f{ri), where the />,• = -iV, (/ = 1, 2, 3) are the generators 
of E(3). Their sum over all i is denoted by K. 

Lemma 1. Define the operator sum T + V by H, with T and V as in 
eqs. (1.1)-(1.2) and assumptions fAlj-fA3j. Let if/ e Ker(E - H) with 
E 6 o'(H), where i// = il/{E;r\,r2,r-i). If V is invariant under the ac- 
tion of E(3), then (a) functions if/ are translation invariant, one writes 
iff = (p(E;ri2,r23,ri2), (b) functions (p solve the eigenvalue equation for 
Tq + V with the same E, namely, 

(To + V')(p(E) = Ecp(E), (2.1a) 

To = -J^(2mir^dl on D(To) (2.1b) 

l<i<3 

where 

5i=Vi2 + Vi3, 52 = V23-Vi2, 53 = -V23-Vi3, 

and(c) D{Tq) = D{T). 

Remark 1. Here and elsewhere, we distinguish potentials V and V by writ- 
ing the equation V,y(r, - rf) = V'.j{rij). Thus V(p = V'ip but V,V = diV . 

Remark 2. Vectors ri2, ^23, ri3 are linearly dependent, ri2 + r23 = ri3, so the 
above given parametrization of (p is rather formal, yet a convenient one for 
our considerations. Hence cp e D(T) c L?-{R^). We shall regard this aspect 
once more in §3. 

Remark 3. It appears that {T + V)if/ = (Tq + V')(p is the identity. For this 
reason, we shall define Tq + V by the same H; there will be no possibility 
of confusion. 

Proof of Lemma L The translation invariance of V infers KV = 0. On the 
other hand, [H,pi\ = -piV. Hence \H,K\ = 0. The three components Kx, 
Ky, of K commute with each other and thus for our purposes, it suflaces 
to choose one of them, say AT^. 

The commutator {H,K^\ = yields if/ G Ker(/l - K^) n Ker(£: - H) 0, 
with A G cr(K;^). Subsequently, functions if/ associated with E are labeled 
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by A as well. One writes ij/ = ipiEX). These functions solve the following 
equation 



The partial differential equation is satisfied whenever functions take 
one of the following forms: exp(i/lzi)^, exp(i/lz2)^ or exp(i/lz3)^, where 
translation invariant functions ip (with K^ip = 0) are labeled by E only, 
(fi = (p{E;r 12, r23,r 12). Since all three forms are equivalent (with their appro- 
priate functions cp), we choose the first one. Note that it suffices to choose 
functions (p being invariant under translations along the z axis only. By ap- 
plying the same procedure for the remaining components, Kx, Ky, we would 
deduce that functions (p are invariant under translations along the x, y axes 
as well. Bearing this in mind, we deduce that functions (p are invariant 
under translations along all three axes associated with each r/y, that is, in 

The application of e(tz) tO(l/ = exp(iAzi)(p yields e(tz)i// = exp(iAtz)il/. This 
means that functions i// labeled by a particular E(3)-scalar representation, 
/I = 0, are translation invariant, il/(EO) = (p(E). 

We wish to find the operator 7^ whose range for all (p in D{Tx) is the 
same as that of T for all in D{T), namely, T;^(p = T\]/. First, we calculate 
the gradients of exp(i/lzi)^. Second, we calculate the corresponding Lapla- 
cians. Third, we substitute obtained expressions in eq. (1.1). The result 
reads 



For a particular /I = 0, we get the tautology Tq = Tq. [Note: di(p = V,(/r.] For 
arbitrary A, substitute in (T + V)il/{EA) = Eifj{EA) and get the equation 
{Tq + V' + AtA)(p(E) = E(p(E). The number of eigenvalues A is infinite, and 
the latter must hold for all of them. It follows from Atx^p{E) = that /l = 
or t^(p(E) = 0. But tx^piE) = is improper since (p = cp(E) is A independent. 
Therefore, A = and functions i// are translation invariant, i/r = (p{E), and they 
satisfy eq. (2.1a). This gives items (a)-(b). Item (c) follows immediately 
due io\]/ = ip. The proof is complete. ■ 

Lemma 2. Define YjiKj'^ij by G, where the sum runs over all I <i < j <3. 
Then there exist domains D, D' c D{Tq) such that D = {(p£ Z)(ro): G(p = 0}, 
DUD' cD(To)andDr)D' = 0. 

Proof. As above, let us choose the component G^. Then the commutators 
[Gz,KJ = 0, [H,KJ = 0, but [G^,H] 0, in general. This indicates that the 



difr dij/ dif/ 
dzi dz2 dz3 
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(only one) eigenvalue A = (proof of Lemma 1) of is degenerate, where 
degenerate eigenfunctions are ^ e D(Tq) (Remark 3). Therefore, if given 
(fi eDi, (pe D', DiUD' = D(To) and Di n D' = 0, then G^ipi = mpi and 
H(f) = Ecp, with some real numbers yu for all ip\ e Di and (peD' . Solutions to 
G^ipi = jjupi are translation invariant functions (pi. In turn, is represented 
by a certain translation invariant function ip\ multiplied by either exp(/izi2) 
or exp(;iZ23) or exp(//zi3). Subsequently, solutions satisfying G^tpx = jjupx 
with // = are in Di as well, by the proof of Lemma 1. The nonempty 
subset D Q Di of these solutions is exactly what we were looking for. ■ 

Remark 4. It follows from the above lemma that H(p = Ecp, with 
H: D(To) L^(R\ actually means H(f) = E(p, with H: D' L^(R^) (see 
also Remarks 2-3). 

Corollary 1. The operator Tq: D ^ I?iR^) is represented by the following 
equivalent forms 



Proof. Combining Lemmas 1-2, simply substitute = in eq. (2.1b) for 



Remark 5. Supposing that all three To-forms are equivalent, we choose the 
first one, eq. (2.2a). Note that, in general. Corollary 1 does not apply to Tq, 
eq. (2.1b), representing the map from D' to L^(R^) (Remark 4). 

2.2. Unitary equivalence. Let Vn, for every integer n = 0,1,..., denote 
G"V', where G" = G^G^ ...G^in times); clearly, is the identity operator 
andVo^ V'- Let/^ (/c = 0,l,...), 

lK = [{rn,r2^,rn)eR'^: rn + r23 = rn,V,+p = 0\/p = l,2,...\, (2.3) 

be a nonempty subset in R^ such that for some arbitrary integer /c > 0, the 
/f* derivative of smooth function V were equal to zero for all k\> k (the 
derivative under consideration is defined in Lemma 2). In what follows. 



Tq=- QrAi2 -;SA23 + 7(^12 • V23), 
= - ^A23 - aAi3 + f(Vi3 • V23), 
= - ^Ai2 -)SAi3 + ?7(Vi2 • V13). 



(2.2a) 
(2.2b) 
(2.2c) 



The real numbers a, J3, 7, ^, ^, rj are equal to 




6eD, and the result follows. 
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we shall identify the spaces endowed with vectors in by k\ hence D^, 
D'l^ and Dq,*;. Here, domains and are considered in a similar way as 
D and D' in Lemma 2 with D(Tq) replaced by Dq ,^, whereas Dq ^ is the 
set of functions from L^il^) such that: (1) ||ro^||^ < oo for all (p e Doy, (2) 
Tq is self-adjoint on Doy, (3) D(Vn) d for all n = 0,1,... Here and 
elsewhere || • ||^ is the L^(/^)-norm. By items (1) to (3), the operator sum 
To + Vn, denoted Hn, is a self-adjoint operator in L^ih) with domain Dq^k- 
In particular, = Tq for all n = /<:+ 1,/f + 2,. . . 

Remark 6. We shall clarify the meaning of L^(/j^)-norm. As it is clear 
from the definition, eq. (2.3), L^(/„) is nothing but L^(R^), with the 
(Lebesgue) measure d/j^ whose exact form depends on the form of V. One 
writes L^(h) = L^{R^,dnid and L^{R^) = L^(R^,diu), where yu^ = //^(/^) and 
c h(R^). Moreover, given k = 0,1,..., the measures {ju^] are mu- 
tually singular since they satisfy jXidR^) = iUk(I k) + /Uk(R^\I k) = fiKih) pro- 
vided ju*:(X\/*;) = by eq. (2.3) for any subset IkCX c R^. Indeed, given 
A: = Oandl. Define /* = /oU/i. Then yUo(/*) +jUi(/*) = //o(^o) +y"i(^*Uo)- 
But /ii(/*\/o) = jJL\{I\). Therefore, for //* = juo+fii, L?-{R^,dix*) is iso- 
morphic to L^iR^ ,diXQ) ® L}{R^ ,diJi), which can be naturally applied to 
arbitrary k. Hence for //' = /io + yUi + . . ., L^(R^,dfi') is isomorphic to 
l}{R^ ,diXQ)®L?-{R^ ,dn\)® . . .. Below we shall demonstrate that fi' = ju. 

Finally, we are in a position to define the L^(/^)-norm. By assumption 
fi' = fi, one may write j^e diJ^ = Y,i^jj d^^, where the sum runs from to oo. 

Equivalently, J/i^ = A,^ Jj^(,d/j. for some nonzero such that Yjk'^k = 1- 

With this definition, that the operator Hn is in L^ih), it actually means that 
the operator A~^Hn is in L^(R^,AKdfi), and therefore there is unitary U such 
that U(A-^Ho@A-^Hi®...)U-'^ is in L^(R^). 

The reason for defining a subset I/^ is to divide a noncommutative re- 
lation [Gz,H] into the commutative ones, as it will be shown thereon. The 
properties of 1^ in the case of Coulomb potentials will be assembled in §3. 1 . 

The main goal of the present paragraph is to demonstrate that the eigen- 
values of H in L^(R^) can be established from the eigenvalues of in 

Theorem 1. Given the operator [u,v]((p) = w((p) in 1^ for k = 0,1,... and 

(f € Dqk. The elements u, v, w denote any operator from G^, Vq, Vi, V^. 
Then('l) 



[G„Vn\ = Vn+i, [Vn,Vm\ = for all n,m = 0,1,2,..., K. (2.4) 
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The element Vn denotes the operator of multiplication by Vn- (2) The com- 
mutation relations in eq. (2.4) define the Lie algebra, denoted ^ = ^(Ik), 
with an operation I^xl^^ I^. 

Proof. First, we shall prove eq. (2.4). Second, we shall demonstrate that 
J?l is a Lie algebra indeed. Note that eventually the commutator [G^, Vn\ 
terminates atn = K due to the definition of I^, eq. (2.3). 

(1) Let us calculate the first commutator in eq. (2.4), namely, 

=(G,V„)(^) = (Gf iy')(^)- 

But (G"'^^V')((p) = Vn+i_((p), by definition. Therefore, the identity is imme- 
diate. The second commutator, [Vn,Vm](^) = VnVm((p) - VmVn(ip) = 0, is 
evident provided Vn and Vm represent numerical functions. 

(2) Elements u - G^, Vq, Vi, ...,Vk form the basis of a linear space Jlih) 
of dimension k + 2. If endowed with the binary operation Jlx Jl ^ Jl, 
denoted (m,v) [u,v] for all u, v in Jl, the linear space Jl must fulfill the 
bilinearity, anticommutativity and Jacobi identity. 

Bilinearity: [au + bv,w] = a[u,w] + b[v,w] for any scalars a, b. Due to 
the commutativity au + bv = bv + au and distributivity au + bu = (a + b)u, it 
suffices to consider two cases: (a) u = Gz,v = Vn,w = Vm', (b) u = G^, v = Vn, 
w = G^; in all cases n,m = 0,l,...,K. 

(a) : [aG, + bVn, VM = (aG, + bVnWm(<p) - V^iaG, + bVn)((p) 

= aGzVm(<p) + bVnVm((p) - aVmGz((p) - bVmV„((p) 
= a[G„Vm\((p) + b[Vn,Vm\((p) 

(b) : [aG, + bVn, GJ(^) = (aG, + bVn)G,((p) - G,(aG, + bVn)((p) 

= aGzGz((p) + bV„Gz((p) - aGfi^{(p) - bG^Vnicp) 
= a[G„G,M<p) + b[Vn,GM 

Hence [•, •] is bilinear. 

Anticommutativity: [m,v] = -[v,u]; in particular, [u,u] = 0. This property 
is easy to verify by using the distributivity of the addition operation. Two 
cases are considered: (a) u = G^, v = V^; (b) a = V„, v = V^; in all cases n, 
m = 0, !,...,/<:. 

(a) : [G„Vnm = G,Vn(<p)-VnG,(<p) = -(VnG,-G,Vn)(<p) = -[Vn,GJ(<p) 

(b) : [Vn,VmM<p) = VnVmW) " V«V„(^) = -{VmVn - V„V«)(^) = -[Vm, V„](^) 

Hence [•,•] is anticommutative. 
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Jacobi identity: [[m,v],w] + [[v,w],m] + [[w,w],v] = 0. The identity is an- 
tisymmetric with respect to the permutation of any two elements. Thus it 
suffices to choose u = G^, v = w = Vm (n, in = 0,1, ... ,k). Then applying 
eq. (2.4) and anticommutativity of [•,•] one finds that 

[[G„VAVM + [[Vn,V^],Gzm + [[V,n,G,],Vnm = [G„Vn]V„^(cp) 

- V^Vn+l((p) + [Vn, V,n]G,((p) - [G„ Vm]Vn(cp) + 
= [Vn+uVM + IVn, Vm+lK^) = 0. 

This completes the proof. ■ 
Corollary 2. The Lie algebra ^{Ik) is nilpotent with the nilpotency class 

K+\. 

Proof. We need only to compute the length of the lower central series con- 
taining ideals = where = Jl, c>l. Assume that c = 1. 
Then Jl^ = is spanned by the elements [m,v], where u, v are in J^. 
Thus by eq. (2.4), the first commutator subalgebra contains elements 
WnYf^^y Similarly, the second subalgebra (c = 2) ^ = is spanned 
by the elements [u,v], where ueJl and v eJl^, hence {VnV^^2 Finally, 
jlK+i _ Q -pj^^ proves that Jl is nilpotent with the central series of length 

K+i. m 

The nilpotency of ^{I,^ implies the existence of an isomorphism from 
S\ to the Lie algebra of strictly upper-triangular matrices ([Hum72, §1.3]; 
[ZS83, §3.1.6]). Equivalently, there exists a representation g: SK^/^) 
given by Q{e) = G^, £)(/„) = V„ (n = 0, 1, . . .,k), where elements g = e,fn e ^ 
denote the strictly upper-triangular matrices; in particular, is isomor- 

phic to the Heisenberg algebra, whereas Jl(Io) is commutative. Clearly, 
[e,fn] = fn+i, [fnjm] = fov all n,m = 0,1, ... ,K, and f„+p = for all p>l. 

Assume now that X is a matrix Lie group with Lie algebra Jl and 
exp: J?l ^ X is the exponential mapping for X. Then the exp(itg) are in 
X for all g e ^ and for all t e R. Provided 11: X ^ GL(7^) is a repre- 
sentation of X on Ii^, we get that n(exp(ig)) = exp(ig(g)). But g(e) = G^, 
and thus n(exp(ie))(0) = I{6), the identity for all 9 e Di^, by Lemma 2. On 
the other hand, g{fn) = V,j (n = 0, 1, . . . ,/c) and n(exp(i/„))(0) = exp(iV„)(0). 
The elements 7, {exp(i?V„)}^^Q therefore form a (bounded) group of uni- 
tary operators given by the map L^ih) for all t e R. In turn, it 
is a subgroup of the group generated by ii7„, where = Tq + y„ and 
Tq: Dq^i, L^ih) is self-adjoint. Due to the filtering J?l D JTI^ D . . . D 
the elements of a set {//n}^=o converge to a single [recall that the ker- 
nel of a Lie algebra homomorphism —> ^'^'^^ is {/c}], which in turn 
commutes with G^, namely, \Gz,Hk^ = 0, by Corollary 2. As a result. 
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the eigenfunctions of operator H^: L^{Ik) those of G^, and thus 

e e Ker(£';, - H,,) n Ker(0 - G^) for e aiH^), by Lemmas 1-2. In 
particular, whenever Iq is nonempty, one should expect that Eq = E due to 
the formal coincidence of Hq with H (Remark 4). However, H is in L^(R^) 
and it is defined on D' whereas Hq in L^{Ik) is defined on at a particular 
value = (see also Remark 5). This means Eq 4^ E, in general (that is, 
for smooth V). On the other hand, provided is nonempty for arbitrary 
K, one finds from the above considered Lie algebra filtering that D' is the 
space decomposition ®kDk, where k goes from to oo. But D' is dense in 
L?-(R^) and in L?-{Ik)- Thus by Remark 6, U is the unitary transformation 



from ®kL\I^) to L\R^) so that UiA'^Hx® A-^H2® . . .)f/~^ = H. Subse- 
quently, for £■ G cr{H) and 6 cr(//^), one finds that £ = Yjk'^kEk, where 
= A;\\\e\\J\mf. But 2.^. = 1 and = ZMl withe = e(E,). One 
thus derives = (||e|U/||</'||)^ and 



As a result, we have established that solutions to the initially admitted 
eigenvalue equation = E^ in L^(R^) are obtained by solving H^d = E^G 
in L^{Ik), where Hk = Tq + Vk with Tq given in eq. (2.2a) and 6 e D^. 

In the next section, we shall be concerned with the Coulomb potentials, 
though one can easily enough apply the method to be presented to other 
spherically symmetric potentials imposed under (A1)-(A3). 

3. Solutions for the three-body Hamiltonian operator with Coulomb 



The Coulomb potential V is a spherically symmetric translation invariant 
function represented as a sum of functions V.j = Zij/rij, eq. (1.2); for the 
notations exploited here, recall Remark 1 . The scalar Zij = Zji = ZiZj, where 
Zi (i = 1,2,3) denotes a nonzero integer (the charge of the particle). The 
spherical symmetry in R^ preserves rotation invariance under SO(3) thus 
simplifying the Laplacian Aij = cfi/drfj + (2 / rij)d / drtj - 1(1 + ^)/ffjy where / 
labels the S0(3)-irreducible representation. Bearing in mind Remarks 2-5, 
we shall use h to label representations for ij =12, and h for ij = 23; the 
associated basis indices will be identified by ;ri = -/i , -/i + 1 , . . . , 0, 1 , . . . , /i 
and by 7T2 = -h, -h + 1, . . . , 0, 1, . . . , /2, respectively. 

3.1. The stability criterion. Let us first study the properties of a subset Ik 
introduced in §2.2. Proceeding from the definition, eq. (2.3), we deduce that 
the nilpotency of is ensured whenever (see also the proof of Corollary 4) 



oo 




(2.5) 



POTENTIALS 
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^ -^=0 for all k = K+p+l=2,3,... for all p=l,2,... (3.1) 

l</<j<3 ij 

provided the z axis is suitably oriented. Equation (3.1) suggests that at 
least one integer Z, from Zi, Z2, Z3 must be of opposite sign - this is what 
we call the stability criterion for the Coulomb three-body system (see also 
Remark 10). There is the classical picture to it: If the three particles are all 
negatively (positively) charged, they move off from each other to infinity 
due to the Coulomb repulsion. A well-known observation follows therefore 
from the requirement that the Lie algebra Jl were nilpotent. Henceforth, we 
accept the criterion validity. 

Linearly dependent vectors {r,j} (Remark 2) form a triangle embedded in 
R^. Based on the present condition, we can prove the following result. 

Lemma 3. Let cj^, o-^, denote the angles between the pairs of vectors 
{^\2' ^u)> (^"13. ^23). (^"12. ^23). respectively. If a given three-body system is 
stable, then for any integer k>2 such that (i) (Z2/Z3) + (Z2/Zi)p^ < and 
(it) 



Ci(^)<- r <— — , Ci{p) = \^ 

(cri is acute) 

I Zi3 \^^^ fO, 0<^<1, 

\ Z12+Z23P''} I—' ^> 1' 

(cTk is obtuse), there exists a multiplier p >0 satisfying sin cr^ = psintok 
so that eq. (3.1) holds for all < Cji>k,cr]^ < n such that: (1) ifO< iOk,crk 
< n/2, then < p < l/sinct»jt/ (2) ifO < a>k< nl2 and nl2 < (Tk <n, then 
1 < ^ < l/sin<t)^; (3) if n/2 < oj^ < n and < (Tk < n/2, then < ^ < 1; 
(4) if n/2 < oJk^cTk < n, then p does not exist. In case that Zi/Zj, < 0, the 
multiplier p {-Z\/Z'i)^l^ for all suitable < 0Jk,crk ^ ^• 

Remark 1. In particular, lemma states that for a certain integer k >2, if 
such exists at all, one can find a multiplier p >0 such that the angles cok, 
CTk, Tk obtained from relations sincTjt = psinojk and sinr/t = Cksinojk (where 
Tk = cok + o-k) solve eq. (3.1). The multiplier Ck = (-Zi^/^Zu +Z23p''))^'''p. 
Clearly, one should bring to mind the sine law relating angles with the as- 
sociated sides of a triangle ri2 + r23 = ri3. 



Corollary 3. Let 
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Dkif) = {0<a)k,(Tk,Tk<n: a)k + o-k-Tk = 0, sin (Tk = p sin ojk, sin Tk = Ck sin ojk]. 

The set Dj^(p) is nonempty ifk>2, the triplet (I, p, c^) fulfills the triangle 
validity, and fi j Zi , Z2 < 0, Z3 > or Zi , Z2 > 0, Z3 < and < -Z\ /Z3 or 
(2) Zi > 0, Z2,Z3 < orZi < 0, Z2,Z3 > and <^ > -Z1/Z3 or (3) Zi,Z3 > 0, 
Z2 < orZi,Z3 < 0, Z2 > and p>0. Otherwise, D^ip) = 0. 

Proof of Lemma 3. Although the proof to be produced fits any positive inte- 
ger k, we shall make a stronger restrictive condition, k>2, due to eq. (3.1). 

The combination of the sine law, r^l sincjjt = ^j"2 ^^^'^k = ^13 sinr^, and 
eq. (3.1) points to the following equation 



[Note that the values ojk, Tk = 0,n are allowed as well by implying crj^ = 0, ;7r.] 
Then the expression for Ck (refer to Remark 7) follows immediately if 
sino-k = psinojk if > 0). The quantity in parantheses 0^^*^ in q is posi- 
tive definite and thus item (i) follows as well. Clearly, the denominator is 
nonzero; otherwise Z1/Z3 < and p (-Zi/Zt,)^^^ must hold. 

By noting that < ri-i/ru = cos£c>^ + (l/p^ - sin^cu^;)^^^, we find from 



where "+" is for 0<crk< n/2, and "-" for n/l <(Tk<n. Items (l)-(4) fol- 
low directly from the above equation: Eg let oJk, (Tk > 7r/2. The denominator 
is of the form (-x)^/*^, x > 0, hence improper (item (4) in lemma). 

Substitute Tk = iOk + cr^ and sin cr^t = psinojk in sinr^: = Qsincy^^ and get 

Ck sin tok = sin(+a;yt + arcsin(g)sina;yt)) 

[here, again, is for <o-k< n/2, and "-" for n/2 <crk<n] yielding the 
estimates 1 - < < 1 + ^ for acute ak, and < Ck < p - I for obtuse ak. 
Provided Ck > 0, substitute the definition for Ck in obtained inequalities and 
get item (ii). This completes the proof. ■ 

Proof of Corollary 3. Items (l)-(3) are obvious due to item (i) of Lemma 3. 
It remains to demonstrate the triangle validity for 1, p, Ck- This is done by 
solving the equation 



^ ^ sin*^ cTyt ^ sin^cr^ ^ 
Z12 + Z23 + Zi3 = 0. 

sm (x)k sm Tk 



eq. (3.1) 



Zl2+Z23S3*^ + 



Zl3 




(cos ojfc + ( 1 / ^2 _ o;^) 1 /2) 
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which yields 

((l+Ck + p)(l+Ck-p)(l-Ck + p)(ck + p-l))^'^ 
sina;^ = (3.2) 

and hence the triangle validity for the triplet (l,p,Ck) must hold due to in- 
equality 0<sincok<l. ■ 

Corollary 4. The spherically symmetric functions can be represented 
by three equivalent forms, where two of them are = Z\2{p;K,k)lr^^ , 
V!{ = Z2-i{p;K,k)lr'^fwith 



Zn(p;K,k) = (-irK\{zn +Z23P'^^ +Zi3(^f ^'), 

andZ23(p;K,k) =Zi2(p;K,k)/p''^^. 

Proof. Differentiate V'^j k times with respect to rij, 

d' Zjj ^ {-IfKlZjj 

Then 

l<;<j<3 ij \<i<j<i ij 

Put into use the sine law and Lemma 3, 

rn_ ^ _ ^ 

^23 ' n3 Ck' 

Substitute the above equations in and get the result. ■ 

Remark 8. Lemma 3 and Corollary 3 provide sufficient information to find 
nonempty subsets 1^. Indeed, consider given nonzero integers Z,- (/ = 1, 2, 3) 
and the real, yet unspecified, multiplier ^ > 0. First, establish possible in- 
tegers k>2,hy Lemma 3. Second, substitute determined values of k in 
Ck- Third, substitute obtained coefficients Ck in eq. (3.2) and get possible 
angles ojk e D^ip) (alternatively, simply apply Corollary 3); the subset 1^ is 
nonempty whenever Di{p) is. Applications to some physical systems will 
be displayed in §3.2.3. 

Remark 9. Although conditions in Lemma 3 and Corollary 3 are invariant 
under the interchange of integers Z, , there might appear some arrangement 
that does not satisfy Lemma 3. If this is the case, one should select another 
one. For example, Zi =Z2 = -1,Z3 = +1, p = 1 brings in Z1/Z3 < and 
p i= I (see Lemma 3), which contradicts the initially defined p = I. On 
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the other hand, Zi = Z3 = -1, Z2 = +1, p = I brings in Ck = 2"^^^, and all 
conditions in Lemma 3 as well as in Corollary 3 are fulfilled. However, if 
none of arrangements of Z, fulfill the lemma, one should conclude that the 
three-body system is unstable. 

Remark 10. We point out that in the present discussion, the definition 
for stability diff"ers from that exploited by [CS90, FB92, Hil77]. Here, 
we do not study the cases of stability against dissociation (see also 
[MRW92, RFGS93]) by assuming these conditions are fulfilled whenever 
bound states are considered. On the other hand, provided the three-body 
system is subjected to the stability criterion, eq. (3.1), one should deduce 
from Lemma 3 that at least Iq 0. As an important example of unbound 
three-body system consider the positron-hydrogen system. Our calculated 
first excited energy [substitute, in atomic units, Zi = Z3 = -l-l, Z2 = -1, 
mi = m2 = 1, m3 = 1836.1527, ni = 2, ^2 = 1 in eq. (3.8b) and then con- 
vert the result into Rydberg units] equals Eq ^ -0.25 Ry, which is almost 
consistent with that of Kar and Ho's [KH05] (see also [DNW78]) derived 
S -wave resonance energy (around -0.257 Ry) associated with the hydrogen 
n = 2 threshold. In the positron-hydrogen system, Ei does not affect the to- 
tal energy E when the first excited states are considered, which means some 
higher eigenstates E^ (/c > 2), if such exist, should be included in order to 
obtain more accurate energies. 

3.2. Eigenstates. We wish to evaluate E e (TdisciH). By eq. (2.5), E is the 
sum of E^ e crdisdHJ, where 

^/. = ^^°+r(Vl2-V23), H^, = Ti+T2 + V' 

and Ti = -aAi2, T2 = -J3A23 (3.3) 

(a, /3 and 7 are as in Corollary 1). We first consider the Hughes-Eckart 
term. Following [Sim70, Appendix 2], we demonstrate that: 

Proposition 1. inf cr(H^) = inf cr(H^). 

Proof. In agreement with Corollary 4 consider in p-space 

= hi, +/3P23 -7(^12 •/'23)' = ap\^ +Zi2(p;K,k)r~^~\ 

[Note that Zi2{f',K,k)r~2~^ can be replaced by Z23{p;K,k)r~^~^; see Corol- 
lary 4. Subsequently, ap^^ is replaced by fip^^ in h,,]. Since /3 > 0, and 
and /SP22 involve independent coordinates, we see inf(r(/?^) = infcr(//0). 

Let p - api2 + bp23 with b > a> Q. Then H',^ = Hk+ jxp^ with /u > 0, 
where //j^ = + (jxab - ylT)q^, where q = pi2 + P23 and 
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H'l^ =hK + \jxa{a -b) + 7/2]pj2 + \Mb(b -a) + y/l^pl^ . 
We choose = yl[la{b - a)] > for & > a > 0. Then fxab - yll equals 
7a/[2(&-a)] >Oand 

=h, + \fi + y{l+bla)ll\ply 

But then mf cr{H'^) = inf cr(/j^.) since and involve independent coor- 
dinates. Subsequently, inf <t{H'i^ = inf cr(H"), for fiab - 7/2 > 0, and finally, 
inf cr{H'i^ = inf <t(H^). Hence inf <t(H^) = inf cr(H^) as desired. ■ 

Remark 1 1. Proposition 1 tells us that the ground state of is that of H^. 

Second, consider H^. It is a self-adjoint operator on Dq^^^ whose eigen- 

functions are in D^. By Remark 6, L^ih) = 1^{R^ ,didK)- But, on the 
other hand, I?{Ik) is isomorphic to I?{R^ ,djiii^\) ® L?'{R^ ^dji^^) provided 
d^K: = dfiK^i ®dfiK,2 [RS80, Theorem 11.10]. We denote L?-(R^,dfiK,i) by 
for / = 1,2. Thus there exist unitary operators Ui and U2 such that 

UiH^U^^ =H^^i®I + I®T2, C/2//^C/2^ = ri(8)7 + /(8)i?j2 
with H^^. = Ti + Vf a =1,2) (3.4) 
and is as in Corollary 4 for / = 1,2. But then 

CTcHsciH^J = (TdisciHl,) = (Td,sc{Hl^) (3.5) 

since cTdisciTi) = crdisciTi) = 0. Equation (3.5) allows one to determine p 
(Lemma 3) for a given k as well as E^^ e <Tdisc{H^)- 

Indeed, an ordinary decomposition of product l}{Q,oo\r^dr)®L}{S^), 
with r = r\2 for / = 1 and r = r23 for / = 2 (5^ is a unit sphere), by an infinite 
sum of S0(3)-irreducible subspaces yields the eigenfunctions O^j of H*^. 

which are of the form C^r~^M^,/(r)y/;r(f2): the normalization constant, 
u^j as in eq. (1.3), Yi„ the spherical harmonics normalized to 1 (recall that 
n = ni for / = 1 and n = n2 for / = 2; the same for / and the spherical angles 
Q). In eq. (1.3), the parameters = Zi2(p;K,k)/a and B^^ = E^/a for / = 1, 
and = Z23{p;K,k)//3 and = £'^/yS for / = 2. Here a and yS (as well as y 
in eq. (3.3)) are as in Corollary 1, and Zi2(p',K,k) and Z23(f,K,k) are as in 
Corollary 4. Therefore, if one solves eq. (1.3) with respect to 5^ for both 
/ = 1 and 2, then the eigenvalues E^ are found from eq. (3.5): E*^ oc B^, 
where the coefficient of proportionality is either a (i = I) or (i = 2). Since 
Bi^ depends on A^ and A^ is a function of p, eq. (3.5) allows one to establish 
p as well. 

Below we shall calculate crdisdH^) and in particular inf (TdisdHK) for in- 
tegers /<■ = and 1. 



ON THE THREE-BODY SCHRODINGER EQUATION WITH DECAYING P0TENTIALS16 

3.2.1. Bound states for k = 0. Allowing Aq < 0, solutions uoj(r) to 
eq. (1.3) appear as a linear combination of the Whittaker [Whi03] func- 
tion VF,j /+i/2(2r ^J-Bo) and its linearly independent, in general, companion 
solution M„^/+i/2(2r V-So), with Bq = -A^I{An^) and n = l+ l,/ + 2,... For 

n>l, Wn,M/2(z) = [(-ir'+H« + 0!/(2/+ l)!]M„,z+i/2(z) and thus M„,,+i/2 
and Wnj+i/2 are linearly dependent. It suffices therefore to select one of 
them, say M„/+i/2(z). The boundary conditions for uoj(r) in L^(0,oo;dr) 
as well as for uoj(r)/r in l?{0,oo\r^dr) are fulfilled: Mnj+i/2iz) and 
Mnj+i/2(z)/z 6io as z ^ 0, and Mn,i+i/2(z) and Mn,i+i/2(z)/z ^ as 



On the other hand, if Ao/r is a repulsive potential, Aq > 0, then uoj(r) 
is represented by a linear combination of functions W_„^/+i/2(2r -^-Bq) and 
M-n,l+i/2(2r V-5o). But M_„,/+i/2(z) ^ oo as z ^ oo and W-n,l+l/2(z)/z 
oo as z ^ 0. Hence none of bound states are observed. However, as pointed 
out by Albeverio et al. [AGHKH04, Theorem 2.1.3] (see also [EFG12]), 
a single bound state exists even if Aq > 0, provided that the Hamiltonian 
operator with / = is defined on a domain of one-parameter self -adjoint 
extensions. 

For an attractive potential, eq. (3.5) yields 



(refer to Lemma 3 for the definition of p) with integers ni = h + lJi+l,... 
and ^2 = ^2 + 1 ^2 + 2, . . . Equation (3.6) indicates that eigenvalues are 
labeled by integers ni , ^2 = 1 , 2, . . . and ^ = 2, 3, . . ., namely, = E{n\, n2,k), 
and 



The procedure to find appropriate k = p + I is described in §3.1 (see Re- 
marks 7-8). In particular, if Docip) is nonempty (Corollary 3), that is, 
for p = oo, eq. (3.1), it holds £■(!,!, oo) < E(ni,n2,oo) < E(n\,n2,k), and 
eq. (3.7) is simplified to 



Z — > oo. 




(3.6) 




(3.7) 
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Remark 12. Assume that given Zi = Z2 = +1, Z3 = -1 and mi < ms. Then 
< f < ni/n2 and k = 2,3,... for all ni < ni- By eq. (3.8a), the lower 
bound of crdisc{H^) equals £(1, l,oo) = -(4jS)"\ which is, under the same 
conditions, in exact agreement with that given by Martin et al. [MRW92, 
eq. (13)] (see also [CS90, §11 A], [FB92, eq. (3a)]). That is to say, for 
the bound three-unit-charge system, inf cTdisdH^) is the lowest bound state 
energy at threshold. 

3.2.2. Bound states for k = 1. We deduce from Corollary 4 that Ai > in 
case Aq < 0. Following the method developed by Nicholson [Nic62], for a 
repulsive potential Ai/r^, we specify the eigenstates which result when this 
potential is cut off by an infinite repulsive core at r = 0. Namely, the only 
solutions uij(r) in L?'(0,oo;dr) which vanish at r = 00 are yfrKy{r ^|-Bl), 

= Al + {I + \ /2)^, where Ky{z) denotes the modified Bessel function of 
the second kind [we specify positive values of v due to Ky{z) = K-y(z)]. 
On the other hand, Mi,/(r) is infinite at r = 0. As demonstrated in [Nic62], 
the solutions yfrKy(r y/Si) exist if r > ro, provided that variables v and 
Bi satisfy Ky(rQ sf-Bi) = 0; ro is known as the cut-off radius. This agrees 
with Case [Cas50] who was the first to establish that for the potentials as 
singular as r~^ or greater, bound states are determined up to the phase factor 
associated with ro. 

The solutions to Mi,/(ro) = are found by expanding Ky(z) in terms 
of Iy{z), the modified Bessel function of the first kind. The result is 
/y(ro yJ-Bi) = I-v(ro V-61) or equivalently. 




(3.8b) 



(3.8a) 



n=l 




00 



[(ro/2)V=gr]^"-^ 
(n-l)\T(n + v) 




Explicitly, 




(r the Gamma function) 



(3.9) 
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for integers n> v, and Bi =0 (no bound states) for integers I <n<v. Sub- 
stitute = Ai + (/ + 1 /2)-^ in n > V and, as in the case when k = 0, come 
by n = / + l,/-i-2,... Subsequently, the coupling constant is bounded by 
<A[ <n^ /2)^ <n^-\/4. It appears that 5i < is unbounded from 

below while the upper bound comes through v ^ n. 

As in the case for k = 0, the eigenvalues e o-disc(H^) are found from 
eq. (3.5), or in particular, from eq. (3.9). The result reads 



Remark 13. The fact that Bi, eq. (3.9), is not bounded from below (for 
n oo) does not necessary mean that cr discing) is unbounded either, as 
this is still to be verified by solving 0-51(1) = /3B[(2) as in eq. (3.10). In 
agreement with Lemma 3, it is apparent that solutions f and k>3 exist only 
for appropriate integers ni and n2 whose range strictly depends on masses 
(or equivalently, on multipliers a, jS). In those cases when none of common 
solutions are obtained, one should deduce that does not affect the total 
energy E and higher eigenstates E^ (k > 2), if such exist at all, should be 
added up to the series of E for obtaining more accurate energies (see also 
Remark 10 for analogous discussion in the case when k = 0). The numerical 
confirmation to it will be given below. 

3.2.3. Some numerical results. To illustrate the application of the approach 

presented in this paper, let us consider the helium atom (He) and the positro- 
nium negative ion (Ps"). Although numerical methods to calculate bound 
states of these physical systems are known in great detail for the most part 
due to Hylleraas [Hyl29], our goal is to comment on results following the 
analytic solutions obtained in the paper. The reason for choosing these 
atomic systems is due to different characteristics of the particles they are 
composed of. We shall calculate some lower bound states associated with 
the scalar SO(3) representation (/ = 0). On that account, the issue of pos- 
sible function antisymmetrization is left out from further consideration as 
well as Proposition 1 holds. All calculations are performed in atomic units. 

The helium atom contains two electrons (Zi = Z2 = -1, mi = m.2 = 1) 
and a nucleus (Z3 = +2, ms = 7294.299536). Here and elsewhere below, 
«i = "2 = 1 (^1 =l2 = 0). First, consider the case /c = 0. Our task is to find 





(3.10) 
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integers k>2 such that Dk{p) 4^ 0. By eq. (3.6), ^ = 0.707155 < 1, thus 
Cfe = ^(1/2 - p^y^^^ exists for all > 3. The variables p, Ck satisfy all nec- 
essary conditions in Corollary 3 for all k>3. Subsequently, vcdadisdHo) 
equals £(1,1, oo) = -2.914048, by eq. (3.8a). Note that the resuk is in- 
variant under the change of charges and corresponding masses (see Re- 
mark 9) with Zi = +1, Z2 = Z-i = -\ and mi = 7294.299536, m2 = ms = 1. 
In this case, p = 1.41412 > 1 and Ck = p(p^/2 - 1)"^^^ exists for all k>3. 
Again, the conditions in Corollary 3 are fulfilled, and the minimal eigen- 
value (refer to eq. (3.8b)) is that obtained just above. In comparison, as- 
sume that Zi = Z3 = -1, Z2 = +2. Then p = I, k>2 and the lower bound 
£(l,l,oo) = -(9/4)q;-i = -4.499383. As seen, the present eigenvalue is 
much lower than the above given one. A somewhat identical tendency is 
observed in all helium- like ions (Li"*", Be^'*', B-^"*" etc.): while the two out of 
three arrangements provide the same eigenvalues, the third one differs and 
it is much lower than the other two. A distinctive feature of this particu- 
lar arrangement is that the multiplier p = I, and vectors {^!j}i</<j<3 form 
the equilateral triangle at ^ = 00, cooo = cToo = tt-Too = 7t/3. To explain 
the appearance of solutions -(9/4)q'~^ ^ -9/2, we refer to the zeroth or- 
der perturbation theory which gives the energy -Z^ = -4, provided that 
the interaction potential r'j'g between two electrons Zi = -1 and Z3 = -1 is 
neglected. This is not the case, as demonstrated above, for the remaining 
two arrangements because both interactions rri and rz} are included in H'^ . 

12 23 K,i 

explicitly; see eq. (3.4). 




(a) The contribution of Ei does not af- (b) The contribution of Ei to the ground 

feet the ground state energy of the he- state energy of the positronium negative 

Hum atom, which is £0 - -2.914048 a.u. ion is well-defined with the inclusion of 

The curves are plotted for k = 3. the cut-off radius. 



Figure 1. (Color online) Solution to aBi(l) = pBi{2) with re- 
spect to k and p; see eq. (3.10). Solutions are found at the points, 
where the plane curves intersect with the same color (the same k) 
dashed curves. 
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Second, consider the case /c = 1. For the arrangement Z\ = Z2 = -1 and 
Z3 = +2, the coefficient = $3(1/2 - ^^Y^'^, hence 1/2 - > for all 
A: > 3. In this case, the estimate 1/2 < vi,y2 < 1 yields A: = 3. However, none 
of < f)< 2"^/^ satisfy ffSKl) = jS5i(2) in eq. (3.10), as it is clear from 
Fig. 1(a). None of common solutions are obtained for the remaining two 
arrangements as well. Following [Cas50], therefore, we deduce that for the 
helium atom, the lowest (ground) eigenvalue of is obtained by the first 
expansion term i?Q in eq. (2.5). 

The positronium negative ion contains two electrons (Z\ = Z3 = - 1 and 
mi = 1713 = 1) and positron (Z2 = +1, m2 = 1); for /c = 0, the other two ar- 
rangements are improper due to Lemma 3: p = (-Zi/Zs)^/*^ = 1 > 2. 
The lowest state inf crdisdHo) is found for ^ = 00 (^ = 1, ojoo = n/S), and it 
equals E(l, l,oo) = -1/4. For /c= 1, the bound 1/2 < vi,V2 < 1 yields p = I, 
k>3. Then vi = V2 = (9/4 -41/-^)!/^ hence k = 3,4,5,6 (see Fig. 1(b)). By 
eq. (3.10), inf o-disciHi) is at = 3 and it equals Ei ^ -0.515488/r5. On 
condition that the ground state of Ps" is -0.261995, by [MRW92], we find 
that the cut-off radius is - 6.56 (eg = 4 in [MP02]). Therefore, for the 
positronium negative ion, only the cut-off" radius ro is needed to calculate 
the ground state energy. 
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